Tetragonal tungsten bronze compounds: relaxor vs mixed ferroelectric - dipole glass 
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We demonstrate that recent experimental data (E. Castel et al J.Phys. Cond. Mat. 21 (2009), 
452201) on tungsten bronze compound (TBC) Ba2Pr2,Ndi_a;FeNb40i5 can be well explained in our 
model predicting a crossover from ferroelectric {x = 0) to orientational (dipole) glass {x = 1), rather 
then relax;or, behavior. We show, that since a "classical" perovskite relaxor like Pb(Mni/3 Nb2/3)03 
is never a ferroelectric, the presence of ferroelectric hysteresis loops in TBC shows that this substance 
actually transits from ferroelectric to orientational glass phase with x growth. To describe the above 
crossover theoretically, we use the simple replica-symmetric solution for disordered Ising model. 

PACS numbers: 77.80.Bh, 77.84.Lf, 77.84.Dy 
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Disordered dielectrics and ferroelectrics for a long time 
have been attracting much attention of scientists due to 
additional (governed by disorder) possibility to control 
their physical properties which might be useful for ap- 
plications. The disordered dielectrics can be in prin- 
ciple divided on two classes. One of them comprises 
the compounds Hke KTaOgiLi (KTL), Nb (KTN), Na, 
which undergo a crossover between ferroelectric and ori- 
entational (dipole) glass phases [H, while the other class 
may be related to so-called relaxor ferroelectrics. Lat- 
ter substances, being actually never ferroelectrics 0,13 j 
belong to perovskite family. Their general formulas are 
Pb(Bi/2B;/2)03 (so-called 1:1 family), Pb(Bi/3B^/3)03 
(so-called 1:2 family) with B-ion being Mn, Zn, Sc, Nb 
or Ta. For example, such materials as Pb(Sc 1/2 Nbi/2 
)03 (PSN), Pb(Sci/2 Tai/2)03 (PST) (1:1 family) or 
Pb(Mni/3 Nb2/3)03(PMN), Pb(Zni/3 Nb2/3)03 (PZN) 
(1:2 family) are considered in the literature to be "clas- 
sical" relaxors. 

The main reason for above relaxors not to become 
ferroelectrics is the nonstoichometry of the position of 
the B ion in their perovskite structure. This nonstoi- 
chiometry actually "spoils" the phonon spectrum inher- 
ent in perovskite structure, destroying the ferroelectri- 
cally important soft mode (see Ref. [j] and references 
therein). Latter fact manifests itself in many observ- 
able quantities of relaxors, to name a few, the smear- 
ing of ferroelectric phase transition and the appearance 
of additional low-temperature peaks in dielectric spectra 
obeying Vogel-Fulcher law There are, however, dif- 
ferences between relaxors and disordered dielectrics like 
KTL. The major difference is that since above relaxors 
do not exhibit macroscopic ferroelectricity, they never 
have the ferroelectric hysteresis, while the substances 
like KTL exhibit it in their ferroelectric and mixed fer- 
roglass phases, see e.g. Ref. [s]. As tungsten bronze 
compound (TBC) Ba2Prj;Ndi_a;FeNb40i5 exhibits fer- 
roelectric hysteresis, which disappears with growth of Pr 
content x, it rather shows the crossower between ferro- 
electric and dipole glass behavior similar to that in KTL 



family. 

In this paper, on the base of analysis of the experi- 
mental facts about TBC 7], we come to conclusion that 
this substance demonstrates mixed ferroelectric- orienta- 
tional glass rather then relaxor behavior. To investigate 
this crossover theoretically, we utilize the simple model, 
based on replica-symmetric solution for disordered Ising 
model. We show that our simple model is able to describe 
qualitatively both the phase diagram and hysteresis loops 
in TBC. Based on our previous results in the KTL fam- 
ily of disordered ferroelectrics, we make predictions about 
dynamical properties of TBC as well as about absence of 
Curie- Weiss law in paraelectric phase of TBC. Latter fact 
may be related to the occurrence of Griffiths (para-glass) 
phase [8] in this substance. The schematic phase dia- 
gram of TBC and similar substances, where a crossover 
from ferroelectric to glassy behavior occurs, is reported 
in Fig. la. This phase diagram is more or less standard 
(see, e.g. Ref. [l|) except the region where Griffiths 
(para - glass) phase is realized. The boundaries of this 
region may, of course, vary from substance to substance 
so that para-glass phase may penetrate deeper in para- 
electric and/or dipole glass phase. This is because the 
different strength of " glassy interaction" between dipoles 
lead to the formation of glassy clusters (which is a physi- 
cal mechanism behind Griffiths phase) at higher or lower 
temperatures. Consequently, these clusters merge into an 
infinite one (signifying the glassy phase onset) at higher 
or lower temperatures. 

To describe theoretically the ferroelectric - orienta- 
tional glass crossover behavior in TBC, we should know 
what kind of dipoles can order ferroelectrically in this ma- 
terial and how many permissible orientations in the crys- 
tal lattice do they have. For example, it is well-known 
(see, e.g. [1] and references therein) that Li in KTa03 
form impurity dipoles due to off-central position of Li in 
a host KTaOs lattice. These dipoles (at T < Tc{x) and 
at a; > Xcr, where a;cr is Li critical content) can order 
ferroelectrically and have six (along [100] kind of direc- 
tions) permissible orientations in a host lattice. As this 
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information is absent in TBC, we describe the thermody- 
namic properties of this substance by the simplest pos- 
sible model of two-orientable dipoles. In such model, to 
account for interplay between ferroelectric and glassy be- 
havior, we use replica formalism for random Ising model. 
Namely, we consider the Hamiltonian 
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where E = Ez is an external electric field (in energy 
units) and the random interactions Jij between pseu- 
dospins Sf and Sj {S^ = ±1) are distributed according 
to Gaussian law 
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where J and A J is, respectively, the mean value (respon- 
sible for long-range ferroelectric order formation) and 
variance (responsible for realization of glassy order pa- 
rameter and destruction of long-range order) of random 
interactions. To describe the dependence of TBC ther- 
modynamic characteristics on Pr content x, the param- 
eters J and AJ should be functions of x. Both these 
quantities will be extracted below from the comparison 
of temperatures Tc (ferroelectric phase transition tem- 
perature) and Tg (glassy freezing transition temperature) 
with their experimental values. 

Conventional replica formalism (see, e.g., [ll,[l3l), be- 
ing applied to Hamiltonian ^ (with respect to Eq. ^j), 
gives following standard (so-called replica-symmetric ^12l . 



13[) equations for long-range m (dimensionless sponta- 
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where ks is Boltzmann constant. These equations (along 
with corresponding free energy function [l3|) define all 
equilibrium thermodynamic properties of the system un- 
der consideration. We note here the well-known fact that 
replica-symmetric solution is unstable against replica- 
symmetry breaking. Under the action of external field E 
one can draw a borderline, separating the regions of sta- 
ble and unstable replica-symmetric solutions. In the vari- 
ables field-temperature, this line is known as de Almeida 
- Touless (AT) line [ij]. Standard (in replica formal- 
ism) but lengthy calculation leads to following criterion 
of replica-symmetric solution stability 



{Ajy 



{kBTfV2^J- 



dy 



cosh 



<1. (4) 



To obtain the AT line from Eq. Q, we should solve it 
(the equality in Eq. (g])) simultaneously with Eqs. ^ for 
order parameters. Our analysis shows that the presence 
of long-range order parameter m plays a role of addi- 
tional (to the external one) stabilizing field so that the 
replica-symmetric solution ([3|) in our case is stable almost 
everywhere on the E — T plane except very low temper- 
ature region. This result is correct for all x except the 
point X = 0, where system is ordered. We note that dif- 
ferent versions of above replica formalism had been used 
to describe the disordered dielectrics like proton glasses 
[lo| and relaxor ferroelectrics (see Ref. [TTl and refer- 
ences therein). 



First we use the equations ([3]) to calculate the T — x 
(T is a temperature) phase diagram of the system, i.e. 
the dependencies Tc{x) and Tg{x), which should be cal- 
culated ioT E = in Eqs. Namely, to find Tc, we put 
TO — )■ in the first equation ^ and expand its integrand 
at small to up to the first nonvanishing term. This yields 
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Substitution of Eq. into the first equation yields 
the vanishing of first term, while the second one gives the 
equation for Tc 
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Now we pay attention to the fact that since g = at 
T > Tg, so that q{Tc > Tg) — 0, we immediately obtain 
from Eq. (O that fc^Tc = J- Then, putting q ^ 
in second Eq. we obtain the expression for Tg{x). 
Finally, the expressions for both transition temperatures 
read 



ksTc = J, ksTq = A J. 



(7) 



To deduce the dependencies of J{x) and AJ{x), we re- 
fer to the Fig. 5 from Ref Q, where the experimental 
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dependencies Tc{x) and Tg{x) for TBC are reported. It 
is seen that there are critical concentrations for glassy 
Xg « 0.2 and ferroelectric a;/ ~ 0.8 orders realization. At 
all other values of x the dependences Tc{x) and Tg{x) are 



Here Jq defines the mean interaction between dipoles 
in TBC at cc = 1, having the physical meaning of a 
mean electric field (in energy units) acting between the 
dipoles. In turn, AJq defines the maximal (at x = 1) 
variance of random interactions. Comparison of Eqs ([7]) 
and ([5]) shows the other meaning of parameters Jq and 
AJq. Namely, Jq defines the mean field (i.e. aX, x < Xg 
where the system is almost completely ordered) ferro- 
electric phase transition temperature, while AJq mea- 
sures the maximal (at x > Xf) glassy freezing tempera- 
ture. The situation here is qualitatively similar to that in 
KTL [Ij. Note that dependences ([5]) reflect the physics 
of the system under consideration: aX x < Xg the width 
of distribution function tends to zero so that Gaussian 
([2]) reduces io 5 - function. It can be shown [3| that 
this 5 - function generates the well-known mean - field 
equation m = tanh(mJo/fcsr) for long-range parame- 
ter m in ordered Ising model. On the other hand, at 
X > Xf the mean value J equals zero so that the dis- 
tribution function ^ remains Gaussian but with zero 
mean, which leads to the situation when the system ad- 
mits only glassy order parameter, see e.g. Ref. [l^. We 
note here that while at all x values except critical con- 
centrations, the dependences Tc{x) and Tg{x) are almost 
linear, near Xf^g the whole power series ([8]) come into 
play. Actually, in Ref. f?\ , the finite values of phase tran- 
sition temperatures appear abruptly near Xf^g. Within 
replica formalism, it is possible to obtain the equations 
for above critical concentrations Xf^g. The calculation of 
Xf^g, however, involves more subtle analysis, then that 
based on simple equations ^ as below Tg the replica- 
symmetric solution ^ is unstable. The results of such 
analysis will be published elsewhere. In the Fig. lb, we 
show the approximation of experimental points (see Fig. 
5 of Ref. f^) Tcix) and Tg{x) by hnear terms in ([8]). The 
fitting linear dependencies were obtained in the form: 

T^x) w 320- 29.049a; (K), 

Tg{x) w 157.2 -f 11.57a; (if). (9) 

From these dependencies we can recover the values Jq ~ 
320/f and AJq ~ 168.8iir. The comparison of these val- 
ues shows that Jq ~ 1.9AJo, which means that mean 
value of interaction is almost two times larger then its dis- 



almost linear. These facts suggest that Pr content depen- 
dence of the phase transition temperatures and hence the 
above parameters of interactions distribution should be 
chosen in the form of power series in x 
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persion. Such large ratio Jq/AJq is actually responsible 
for the situation when TBC exhibits ferroelectric long- 
range order for the majority of Pr contents - from a; = 
to a; = 0.8. It is seen that the above linear dependencies 
approximate the experimental points pretty good. This 
means that our simple replica-symmetric solution cap- 
tures qualitatively all peculiarities of a crossover between 
orientational glass and ferroelectricity in TBC. 

Next we use the equations ^ for the qualitative de- 
scription of ferroelectric hysteresis m{E) in TBC. Fig. 2a 
reports the qualitative picture of a hysteresis, described 
by the equations ([3]). Namely, if T > Tc{x) (paraelectric 
phase), the curve m(E) is monotonously growing func- 
tion so that both directions of electric field variation are 
the same, i.e. we are on the same curve. At T < Tc{x) 
the m{E) curve acquired "s-shape" with central unsta- 
ble part with dm/dE < 0. This unstable part yields the 
hysteretic behavior, namely jumps at _B = iEc, where 
Ec is a coercive field. These jumps, as usual, generate 
the hysteresis loop as shown on Fig. 2a. The results of 
specific calculations for the experimental tempera- 
ture T = 280K (corresponding to dimensionless quan- 
tity ksT/Jo = 0.875) are reported in Fig. 2b. It is seen 
that maximal coercive field is achieved at a; = 0, where 
the whole system is ferroelectric. As Pr content x grows, 
the coercive field decays so that at x corresponding to Tc 
value, the unstable part of the m{E) curve degenerates in 
a vertical line and Ec — 0. At x values corresponding to 
T > Tc, we have the paraelectric phase with monotonous 
m{E) curve similar to that from Fig. 2a. We note that 
saturation polarization diminishes as x grows. This is 
because at larger x the dilution of the system grows so 
that there is less (then at a; = 0) ferroelectrically active 
dipoles, which give smaller saturation polarization. 

The comparison of our theoretical hysteresis loops from 
Fig. 2 with experimental ones from Ref. 0] shows their 
different shapes in theory and experiment. While both 
theory and experiment give the decay of saturation po- 
larization as X grows, the coercive field in the experi- 
ment and theory have different behaviors. Namely, the 
experimental coercive field depends on x only weakly, 
while the theoretical one strongly. Latter fact gives the 
different shapes of theoretical and experimental hystere- 
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FIG. 1: The phase diagram of the system under consid- 
eration. Panel (a) depicts schematic picture with Griffiths 
phase shown. Panel (b) reports the approximations Tc(x) — 
290.95 + 29.05(1 -a;) KandTgix)^ 157.2 + 11.57a; (full lines) 
and experimental dependences [3] Tc{x) and Tg{x) (squares). 
Critical x values for glassy Xg and ferroelectric Xf are also 
shown. 



FIG. 2: The theoretical hysteresis loops m{E/Jo), obtained 
from solution of equations of Eqs. ((3]). Panel (a) shows quali- 
tative situation. Dashed line corresponds to the unstable part 
of the m{E/ Jo) curve at T < Tc{x). Arrows show the possible 
directions of sweep around hysteresis loop. The coercive field 
Ec is also shown. Panel (b) reports quantitative solution for 
ksT — 0.875Jo and different x. 



sis loops. To calculate the above shape more precisely, 
the consideration of repolarization of domain structure 
of TBC in the external electric field is necessary. This 
problem requires experimental investigations of character 
of ferroelectric domain structure in TBC. The theoreti- 
cal approach dealing with domain structure characteris- 
tics in disordered ferroelectrics, have been put forward 
earlier 16|. Other way of (less physical) calculation of 
hysteresis loops in TBC is to use Preisach model [l^l - 
like approach (see e.g. [l8| and references therein), where 
the loops m{E) from Fig. 2b (we recollect that these are 
actually the solutions of Eqs. ^ at T < Tc{x)) can be 
considered as elementary hysterons. As those hysterons 
are functions of T and x, the resulting loop would also 
depend on these parameters. 

To conclude, in the present paper we have 
shown that tungsten bronze compound (TBC) 



Ba2Pr2;Ndi_a:FeNb40i5 undergoes a crossover be- 
tween ferroelectric and orientational (dipole) glass 
phases rather then exhibits "classical" (i.e. that adopted 
in the literature) relaxor behaviour. To demonstrate 
that, we use the replica-symmetric solution of disordered 
Ising model which permits to derive the dependences of 
phase transition temperatures Tc and Tg on Pr content x. 
We note here that the above replica formalism permits 
to obtain the equations for critical concentrations for 
ferroelectric Xf and glassy order Xg appearance. This 
issue can be elaborated considering the nonequilibrium 
(like frequency dependent dielectric susceptibility) 
properties of the system under consideration; it is out of 
frames of present paper. 

The same replica-symmetric equations ^ have been 
used to calculate the dependence m{E), determining 
qualitatively the ferroelectric hysteresis in the system. 
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We have shown that hysteresis behaves according to the 
system phase diagram. Namely, as (at given x) T ap- 
proaches Tc or a; — > x/, the loops monotonously shrink, 
giving smaller coercive field values. Ki T ~ Tc ov x — x f 
the loop has zero width and coercive field Ec = 0- Al- 
though our theory describes the decay of saturation po- 
larization at Pr content growth, we were not able to ex- 
plain the enhancement of Ec (the "width" of hysteresis 
loop) in TBC at a: = 0.6, see Fig. 6 of Ref. 0. The ex- 
planation of this puzzling behavior requires (even on the 
level of replica-symmetric solution ([3])) more exquisite 
approaches (then above simple model) like supposition 
that there are additional defects (interacting with fer- 
roelectrically active dipoles) in the system such that at 
certain x (e.g. &i x — 0.6 in Ref. 01) they enhance 
the spontaneous polarization and consequently coercive 
field. The experimental fact |J7|] that paraelectric phase 
of TBC does not follow the Curie- Weiss behavior, can be 
explained in terms of Griffiths (paraglass) phase realiza- 
tion in the substance Q, see Fig. la. Although the dy- 
namics of order parameters cannot be calculated with the 
help of above replica-symmetric solution, its qualitative 
description can be made similar to that in ferroelectric 
polymers [13] ■ Such calculations are out of frames of a 
present work and will be published elsewhere. We em- 
phasize once more that to elucidate the interplay between 
ferroelectric and orientational glass behavior in nonper- 
ovskite disordered ferroelectrics, additional experimental 
investigations are highly desirable. 

This work was supported in part by Opole University 
Intramural Grant (Badania Statutowe). 



[1 
[2 

[3; 

[4: 

[5: 
[6: 
[t: 
[s: 

[9: 
[lo: 

[11 

[12 

[13 

[14 

[15 
[16 

[17 

[is: 

[19 



U.T. Hochli, K. Knorr and A. Loidl Adv. in Physics 39, 
405 (1990). 

G.A. Smolenskii, A.I. Agranovskaya Sov. Phys. Sol. State 

I, 1429 (1959). 

W. Dmowski, S. B. Vakhrushev, I.-K. Jeong, M. P. 
Hehlen, F. Trouw, and T. Egami Phys. Rev. Lett. 100, 
137602 (2008). 

P. M. Gehring, H. Hiraka, C. Stock, S.-H. Lee, W. Chen, 
Z.-G. Ye, S. B. Vakhrushev, and Z. Chowdhuri, Phys. 
Rev. B 79, 224109 (2009). 

D. Vieland, S.J. Jang, L.E. Cross and M. Wuttig, J. Appl. 
Phys. 68, 2916 (1990). 

L.A. Boatner, U.T. Hochli, H. Weibel, Helv. Phys. Acta 
50, 620 (1977). 

E. Castel, M. Josse, D. Michau and M. Maglione J. Phys. 
Cond. Mat. 21, 452201 (2009). 

V.A. Stephanovich European Phys. Journ. B, 18, 17 
(2000). 

B. E. Vugmeister and M. D. Glinchuk, Rev. Mod. Phys. 
62, 993 (1990). 

R. Pirc, R. Blinc, V. Bobnar, and A. Gregorovic Phys. 
Rev. B, 72, 014202 (2005). 

R. Pirc, B. Tadic, R. Blinc, Phys. Rev. B, 36, 8607 
(1987). 

LYa. Korenblit, E.F. Shender Sov. Phys. Uspekhi, 32, 
139 (1989). 

K. H. Fischer and J. A. Hertz, Spin Glasses, (Cambridge 

University Press,. Cambridge, 1991). 

J.R.L. de Almeida and D.J. Thouless, Journ. Phys. A, 

II, 983 (1978). 

V.A. Stephanovich Ferroelectrics, 192, 29 (1997). 

E. V. Kirichenko and V. A. Stephanovich, Ferroelectrics, 

390, 99 (2009). 

P. Preisach, Zeitschrift fiir Physik, 94, 277 (1935). 
E. Delia Torre, Physica B, 343, 1 (2004). 
V.A. Stephanovich, M.D. Ghnchuk, E.V. Kirichenko, B. 
Hilczer, Journ. Appl. Phys. 94, 5937 (2003). 



* URL: http : //draco .uni . opole .pl/~ stef an /VStephanovichPossier .htmi] 

Electronic address: stef@math. uni. opole. pl^ 



